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We study theoretically a radio frequency superconducting interference device integrated with a 
nanomechanical resonator and an LC resonator. By applying adiabatic and rotating wave approxi- 
mations, we obtain an effective Hamiltonian that governs the dynamics of the mechanical and LC 
resonators. Nonlinear terms in this Hamiltonian can be exploited for performing a quantum nonde- 
molition measurement of Fock states of the nanomechanical resonator. We address the feasibility of 
experimental implementation and show that the nonlinear coupling can be made sufficiently strong 
to allow the detection of discrete mechanical Fock states. 

PACS numbers: 03.65.Yz, 85.25.Dq 



Anharmonicity may introduce coupling between differ- 
ent modes of a resonator. In a seminal paper Sanders 
and Milburn found that inter-mode coupling could in 
principle enable quantum non-demolition (QND) detec- 
tion 0, H| of discrete Fock states of a signal mode by 
intensively driving another nonlinearly coupled detector 
(or pump) mode, and monitoring the response near the 
pump frequency. Such a measurement scheme is charac- 
terized by a measurement time r m , defined as the time 
needed to distinguish between initial states of the signal 
mode having different Fock numbers. Fock number de- 
tection can be realized if r m can be made shorter than 
the lifetime of a Fock state of the signal mode. 

The prospects of employing this technique for QND 
measurement of a single phonon in a mesoscopic mechan- 
ical resonator was recently studied theoretically in Ref. 
0. In practice, however, anharmonic couplings between 
different modes of a mechanical resonator are typically 
far too weak to allow single-phonon detection when the 
detector mode is taken to have a linear response. On the 
other hand, a significant enhancement can be achieved by 
driving the detector mode into the nonlinear regime and 
exploiting critical slowing down by operating close to the 
edge of the region where the response exhibits bistability 
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FIG. 1: The device consists of an RF SQUID integrated with 
a nanomechanical and an LC resonators. The external flux is 
applied by using 2 bias current lines carrying dc currents hi 
and 7i,2 respectively. One of the current line lies very close 
to the vibrating beam, allowing thus a relatively high value 
of the magnetic field B at the location of the beam. On the 
other hand, the Josephson junction is located symmetrically 
in such a way that the magnetic field at its location vanishes 
when hi = It2- This configuration allows applying the needed 
external flux to the loop <I>e = $o/2 without significantly 
degrading the critical current of the junction provided that 
the area of the junction is kept much smaller than the area of 
the loop. 



To achieve single phonon sensitivity, it is highly desir- 
able to increase the anharmonic coupling between the sig- 
nal mode and the detector mode. Here we propose a novel 
configuration in which a fundamental mode of a doubly 
clamped beam, serving as the mechanical signal mode, 
is indirectly coupled to an LC resonator serving as the 
detector. The anharmonic coupling is achieved by cou- 
pling both the mechanical mode and the LC resonator di- 
rectly to a radio frequency superconducting interference 
device (RF SQUID) Q. We show that the RF SQUID 
degree of freedom can be eliminated from the equations of 
motion by employing an adiabatic approximation. More- 
over, employing the rotating wave approximation (RWA) 
leads to further simplification of the equations of motion 



of the system. We find that the effective anharmonic cou- 
pling constant between the mechanical mode and the LC 
resonator in this configuration can be made sufficiently 
large to allow single phonon detection under appropriate 
conditions. 

Alternative approaches for performing QND measure- 
ment of the energy of a nanomechanical resonator were 
proposed recently in Refs. [8|, [9j. Similar systems con- 
sisting of a SQUID integrated with a nanomechanical res- 
onator have been recently studied theoretically [13, EH • 
Zhou and Mizel have shown that nonlinear coupling be- 
tween a DC SQUID and a mechanical resonator can be 
employed for producing squeezed states of the mechanical 
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resonator [12| . More recently, Xue et al. have shown that 
a flux qubit integrated with a nanomcchanical resonator 
can form a cavity quantum electrodynamics system in 
the strong coupling region (l3j . 

The device, which is seen in Fig. [TJ consists of an 
RF SQUID inductively coupled to an LC resonator. A 
section of the RF SQUID loop is freely suspended and 
allowed to oscillate mechanically. We assume the case 
where the fundamental mechanical mode vibrates in the 
plane of the loop and denote the amplitude of this flex- 
ural mode as x. Let to be the effective mass of the fun- 
damental mode, and Lu m its angular resonance frequency. 
A magnetic field is applied perpendicularly to the plane 
of the loop. Let $ e be the externally applied flux for 
the case x = 0, and let B denote the component of the 
magnetic field normal to the plane of the loop at the lo- 
cation of the doubly clamped beam (it is assumed that 
B is constant in the region where the beam oscillates). 
A Josephson junction (JJ) having a critical current I c 
and capacitance Cj is integrated into the loop. The RF 
SQUID is inductively coupled to a resonator comprising 
an inductor L and a capacitor C in parallel. The mutual 
inductance between the RF SQUID and the resonator is 
M. Detection is performed by injecting a monochromatic 
input current ij n into the LC resonator and measuring 
the output voltage V out (see Fig. [TJ). 

The total magnetic flux <!> threading the loop is given 
by $ = $ e + Blx + $i, where I is an effective length 
of the beam. The term $i = I S A + MIl represents the 
flux generated by both the circulating current in the RF 
SQUID I s and by the current in the inductor of the LC 
resonator II, where A is the self inductance of the loop. 
Similarly, the magnetic flux in the inductor of the LC 
resonator is given by <p = II L + MI S . 

The gauge invariant phase across the Josephson junc- 
tion is given by 6 — 2irn — 271"$/ $o, where n is integer 
and $o = h/2e is the flux quantum. 

The Lagrangian of the closed system is expressed as a 
function of x, (p and <& and their time derivatives (denoted 
by overdot) 



The Euler - Lagrange equations can be written as 



„ mx 2 Cip 2 Cj$ 2 TT TT 
where the potential terms are given by 
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(3a) 
(3b) 

(3c) 



The interpretation of these equations is straightforward. 
Eq. (|3a[) is Newton's 2nd law for the mechanical res- 
onator, where the force is composed of the restoring elas- 
tic force —mco^x and the Lorentz force Bll s acting on 
the movable beam. Eq. (]3b[) states that the injected cur- 
rent Ii n into the LC resonator equals the sum of the cur- 
rent in the inductor II and the one in the capacitor C(p. 
Similarly, Eq. (|3c[) states that the circulating current I s 
equals the sum of the current I c sin 8 through the J J and 
the current CjVj through the capacitor, where the volt- 
age Vj across the JJ is given by the second Josephson 
equation Vj = ($o/27r)0. 

The variables canonically conjugate to x, ip and $ are 
given by p — dC/dx, q = d£/dip and Q = d£/d& re- 
spectively. The Hamiltonian is given by 



where 



Ho = 
Hi = 
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(5a) 
(5b) 



Quantization is achieved by regarding the variables x, p, 
(p, q, 4> and Q as Hermitian operators satisfying canonical 
commutation relations. 

As a basis for expanding the general solution we use 
the eigenvectors of the following Schrodinger equation 



Hi \n (x, <p)) = e„ (x, ip) \n (x, ip)} 



(6) 



where x and ip are treated here as parameters (rather 
than degrees of freedom). The local eigen- vectors are 
assumed to be orthonormal (to (x, ip) \n (x, ip)} = S nm . 

The eigenenergies e n (x, ip) and the associated wave- 
functions #„ are found by solving the following 
Schrodinger equation 
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where 



muj 2 n x 2 CJjip 2 
Uq = — o 1 hri'P , 
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and where uj e — 1/ \JLC and K = Mj \J~KL. 
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L = 27rA/ c /$ , (3c = 2e 2 /CjE and E Q = 
The total wave function is expanded as 



*§/8tt 2 A. 



£ n (x,ip,t) \n) 



In the adiabatic approximation [14j | the time evolution 
of the coefficients is governed by the following set of 
decoupled equations of motion 



[Ho 



(x, ip)] £„ = ih£r, 



(11) 



Note that in the present case the geometrical vector po- 
tential [l4| vanishes. The validity of the adiabatic ap- 
proximation will be discussed below. 

In what follows we focus on the case where \(f>o\ <C 1 
and 0l (l — K 2 } > 1. In this case the adiabatic poten- 
tial u (</>) given by Eq. (JSJ contains two wells separated 
by a barrier near <j> = 0. At low temperatures only the 
two lowest energy levels contribute. In this limit the lo- 
cal Hamiltonian Hi can be expressed in the basis of the 
states \-r\) and \rx), representing localized states in the 
left and right well respectively having opposite circulat- 
ing currents. In this basis, Hi is represented by the 2x2 
matrix 



V<t>o A 

A -T](j) 



(12) 



The real parameters 77 and A can be determined by 
solving numerically the Schrodinger equation {?]). The 
eigenenergies are given by 



(13) 



It is convenient to introduce the annihilation operators 



A m — 



K. = 



y/muj m x - 



V Cu> e ip + 



y/CuJe 



(14) 
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and the corresponding number operators N m — A m A m 
and N e = A\A e . Consider the case where $ e = $o/2, 
and assume that adiabaticity holds and the RF SQUID 
remains in its lowest energy state. The energy of this 
state £_ (2, ip) can be expressed in terms of the annihi- 
lation operators A m , A e and their Hermitian conjugates. 
The resulting expression generally contains terms oscil- 
lating at frequencies u>k,i = kuo m + lu> e , where k and I are 
integers. In the RWA such terms are neglected unless 
Uk.i — or u>k,i is close to the frequency of the externally 
injected bias current Ii„, since otherwise the effect on the 
dynamics on a time scale much longer than a typical os- 
cillation period is negligibly small. Keeping terms up to 
fourth order in <f)Q the eigenenergy e_ in the RWA can 
be expressed as 



£rwa (a;, <p) = H^o + &-2,oN m + tt , 2 N e 

+ n^ N m + q qa n 2 + n 2a N m N e ) 



(16) 



The constant term Mlo can be disregarded, since it only 
gives rise to a constant phase factor. Moreover, the lin- 
ear terms Ml 2 ,oN m and h£lo. 2 N e only give rise to a small 
renormalization of the resonance frequencies uj m and Lo e 
respectively. In addition, since no external drive is ap- 
plied directly to the mechanical oscillator, the expecta- 
tion value (N m ) is expected to be relatively small and 
consequently the term proportional to N 2 n can be ne- 
glected. Thus, in the RWA the Hamiltonian Ho + e_, 
which governs the dynamics of the mechanical and LC 
resonators, is given by 

H.rwa = hu m N m + hu e N e 



2Cu e K 

+ nn 0A N 2 + nn 2 , 2 N m N e , 



where 
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(18) 
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(19) 

All subsequent analysis is based on Refs. [H, 0, S @], 
which consider how nonlinear terms of the kind appear- 
ing in the Hamiltonian 7Yrwa can be utilized for QND 
detection of discrete Fock states. A single phonon added 
to the mechanical resonator results in a shift in the effec- 
tive resonance frequency of the LC resonator. Consider 
first the case f2o,4 = 0, where the response of the LC res- 
onator is linear, and thus this frequency shift is given by 
^2,2- The effective resonance frequency can be contin- 
uously monitored by employing an homodyne detection 
scheme, in which V ou t is mixed with a local oscillator at 
the same frequency as the frequency of the driving cur- 
rent Ii n . The measurement time r m , which is needed to 
detect a frequency shift ^2,2, is given in the high tem- 
perature limit k B T > hw e by P. Iiol [TH 



27T7e ksT 



^2,2 



(20) 



where Uq — (N e ) fku e is the stored energy in the stripline 
resonator. Single phonon sensitivity is achieved when 
£ > 1, where the dimensionless parameter £ is defined 
as £ = to/T m , and to is a characteristic lifetime of a 
Fock state of the mechanical resonator. At a high tem- 
perature ksT ^> hw m the lifetime is given by tg 1 ~ 
7 m (ksT '/hoj m ) 2 , where j m is the damping rate of the 
mechanical mode and ksT is the thermal energy (see 
Eq. (55) of Ref. ||). 

When, however, Sl ,4 7^ 0, the response of the LC res- 
onator is approximately linear only when (N e ) is smaller 
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than the critical value corresponding to the onset of 
nonlinear bistability. Using this critical value (N e ) — 
7 e /v / 3^o,4, which is given by Eq. (38) of Ref. one 
finds that the largest possible value of C in the linear 
regime is roughly given by 



_ 12A / j] 2-kBI / h \ fhuj m \ hui e 



2.8 x 10~ 15 x -rA-^ f — 



h GHz 7,„ V A T (im j 2tt GHz 



2ttGHz 



( 10-" kg ) (k) 



(21) 



As was mentioned above, operating the LC resonator 
pump mode in the regime of nonlinear response may al- 
low a significant enhancement in the sensitivity by driv- 
ing the LC resonator to a jump point at the edge of the 
region of bistability 0, i, [lCj. Note, however, that the 
perturbative approach employed in Ref. [5] is not valid 
close to the edge of the region of bistability, and thus 
further study is required to analyze the behavior of the 
system in this region. 

We now return to the adiabatic approximation and ex- 
amine its validity. As before, consider the case where the 
externally applied flux is given by <& e = $o/2, and the 
LC resonator is driven close to the onset of nonlinear 
bistability where the number of photons approaches the 
critical value (N e ) . Using Refs. [7J, ll6| one finds that 
adiabaticity holds, namely Zener transitions between adi- 
abatic states are unlikely, provided that 



ttA 2 



> 1 



where 



2nM 2(N e ) c hoj e 



C 



(22) 



(23) 



As can be seen by examining Eq. (|2"Tj) , satisfying the 
condition ( max ^ 1 together with all other requirements 
is quite difficult. However, as we show below, a care- 
ful design together with taking full advantage of recent 
technological progress may open the way towards experi- 
mental implementation. We discuss below some key con- 
siderations. Due to the relatively strong dependence of 
(max on r;/A and on Bl, it is important to minimize the 
uncertainty in the values of these parameters to allow 
a proper design. The parameter Bl plays a crucial role 
in determining the coupling strength between the me- 
chanical resonator and the RF SQUID. Enhancing the 
coupling can be achieved by increasing the applied mag- 
netic field at the location of the mechanical resonator B. 
However, B should not exceed the superconducting criti- 
cal field. Moreover, the externally applied magnetic field 
at the location of the JJ must be kept at a much lower 
value in order to minimize an undesirable reduction in 



I c . This can be achieved by employing an appropriate 
design in which the applied field is strongly nonuniform 
(see Fig. [T]). In addition the sensitivity of the device can 
be further enhanced by increasing the ratio rj/A. To in- 
crease the value of r\ it is desirable to employ a JJ having 
a high plasma frequency. On the other hand, the value of 
the energy gap A can be reduced by increasing [3l- Note, 
however, that A has to be kept larger than fc^T to ensure 
that thermal population of the first excited state of the 
RF SQUID is negligible. In addition, a successful exper- 
imental implementation requires a very low temperature 
T, low mass m, high resonance frequency cj e , and low 
mechanical damping 7 m . Note also that, as was pointed 
out before, an additional enhancement of ( (beyond the 
linear value (max) may be achieved by exploiting nonlin- 
earity 0, la, LlOl • 

As an example, we consider below a device in which 
a single-walled carbon nanotube, which serves as the 
nanomechanical resonator UtI. 18l. H^l . is integrated into 

To en- 
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the loop of an RF SQUID 
hance the plasma frequency we assume the case where a 
microbridge weak link is employed as a Josephson junc- 
tion in the RF SQUID [25J]. The following parameters 
are assumed: C = 3.2 x 1CT 11 F, w e /2vr = 0.22 GHz and 
uj e / le = 10 4 for the LC resonator, A = 1.1 x 10~ 10 H, 
Cj = 10- 16 F and I c = 12 fiA for the RF SQUID, 
m = 10~ 19 kg, iv m /2n = 0.5 GHz and uj m /j m — 10 3 for 
the mechanical resonator, temperature T — 0.02 K, and 
K = 0.001 and Bl = 0.05 Tx /im for the coupling be- 
tween the RF SQUID and between the LC and mechan- 
ical resonators respectively. These example parameters 
are achievable with present day technology. The chosen 
value of A corresponds to a circular loop with a radius 
of about 18 /im and a wire having a cross section radius 
of about 0.22 fim, whereas the values of Cj and I c corre- 
spond to a junction having a plasma frequency of about 
1.2 THz. Using these values one finds (3l = 3.9 and (3c = 
1.0. The values of (3l, Pc and K are employed for cal- 
culating numerically the eigenstates of Eq. ([7]) [7f • From 
these results one finds for the values of the 77 and A pa- 
rameters in the two-level approximation to Hamiltonian 
Hi [Eq. (TT3])] 77 = 4.8^0 and A = 1.2 x lO" 3 ^- Using 
these values yields ^2,2/2^ = 220 kHz, A/h = 0.92 GHz, 



Cmax — 230 

ttA 2 

170 



rjhT c 



(24a) 
(24b) 



Eq. (I24ap indicates that single phonon detection is fea- 
sible, even when the response of the LC resonator is as- 
sumed to be linear. Moreover, Eq. (|24b[) ensures the 
validity of the adiabatic approximation. 

In summary, we have demonstrated that an RF SQUID 
can be exploited for introducing a tailored coupling be- 
tween a mechanical and LC resonator, and that such 
coupling can be made sufficiently strong to allow a QND 
measurement of discrete mechanical Fock states. Accord- 
ing to the Bohr's complementarity principle [26| . in the 
limit of single-phonon sensitivity, dephasing is expected 
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to come into play, leading to broadening of the resonance 
line shape of the mechanical signal mode, allowing thus a 
relatively simple experimental verification by performing 
a which-path like experiment. 
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